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RHIC experiments have established jet quenching as a fundamental tool in the study of hot 
' matter in heavy-ion collisions. The energy reached is enough for high-pT inclusive particle 

studies and one- or two-particle correlations up to px — W — 20 GeV. The corresponding data 
^ • are successfully described by a formalism based on radiative energy losses. This formalism 

^ ' needs, however, to be improved for more exclusive observables, in particular for the recon- 

,.1^ , structed jet properties. Jet reconstruction is a new branch in heavy-ion collisions with very 

high expectations as a tool for the characterization of the medium properties. We present 
here an implementation of the jet quenching phenomena based on a modification at the level 
. of the splitting functions in parton showers in standard, general purpose, Monte Carlos as 

■ PYTHIA or HERWIG. 

in 

C5 I Jet physics is a recent subject in heavy-ion colHsions. Indeed, RHIC (and partially also 

' the SPS) have shown in the last ten years how rich the physics studies at the largest possible 

■ transverse momentum of the produced particles can be^H. Most of those data correspond to 
, inclusive one- or two-particle measurements, where trigger bias effects appear. These effects can 

be generically understood as a geometric bias in which only those particles close to the surface 
of the produced medium can escape it. Reconstructed jets have long been regarded as ideal 
tools to surmount this difficulty: in a perfect jet reconstruction no energy is "lost" and the bias 
d I associated to triggering in a steeply-falling px spectrum disappears. The main difficulty in such 

studies is to deal with the large multiplicity environment of a heavy-ion collision in which the 
underlying event represents ~ 100 GeV per unit area at RHIC — and ~ 250 GeV is expected 
at the LHC — i.e. figures similar to the total energy of the eventually reconstructed jets. This 
difficult analyses benefitted, ho wev er, from tools initially developed to control de pileup in the 
LHC proton-proton programme'^lSI which found a natural application in heavy- ion collisions. In 
this way, reconstructed jets at RHIC were first presented at the Hard Probes conference in 2008 
1^1^ and, latter on and with important improvements, at Quark Matter 2009 and subsequent 
conferences. Full control over different experimental biasses and procedures is under ongoing 
studies by the experimental collaborations. These works represent the beginning of a completely 
new avenue in the study of heavy-ion collisions. 

On the theory side, the description of experimental data in terms of the so-called radiative 
energy loss is very successful — see e.g. Ref. ^ for recent reports on the subject — although 
some issues are still open, the most important one being the suppression of heavy quarks in 
the hot medium. All phenomenological implementations rely, in one way or the other, on 
the independent gluon emission approximation, first proposed in Ref. ^ precisely to deal with 
the bias effects described above. This prescription is enough to describe inclusive one- or two- 



X 



particle production due to the fact that the medium-induced radiative spectrum is not divergent. 
However, the way to describe more exclusive observables as jet structures or intra-jet correlations 
is not known from first principles calculations. Several attemp ts exist to deal with this problem, 
in most cases by Monte Carlo implementationl- ^ l ^ l -*'^ l -*'-*^ l -*^2 | 13 | 

We present here a modification of the PYTHIA parton shower routine to include the 
final-state medium-modified parton shower evolutional^] similar modification in 
is under way^^. 



1 A medium-modified parton shower evolution 

Gluon radiation, independent except for an ordering variable (virtuality, angle, etc.), leads to 
an excellent description of jet evolution in the vacuum. The ordering is dictated by the presence 
of divergencies in the splitting probabilities which are resummed in well-controlled approaches. 
Monte Carlo generators such as PYTHIA or HERWIG encode this physics in parton shower 
routines. So, gluon multiplication is the main building block of jet structure formation in the 
vacuum. 

Several different effects can affect the jet evolution in the presence of a medium, among 
others: 



1. Medium-induced gluon radiation — which modifies the splitting probability producing 
broadening, softening of the spectrum and enhancement of the jet multiplicity. 

2. Non-eikonal corrections to the splitting probability — collisional energy loss producing a 
flow of energy (dominantly) from the fast partons to the medium. 

3. Modification of the color structure of the jet evolution — exchanging color with the medium 
would lead to a different color structure at the end of the parton shower, modifying the 
hadronization stage. 

4. Role of the ordering variable — the space-time picture of the jet is irrelevant in the vacuum 
but not in the medium where an interplay between the extension of the medium and that 
of the jet could be present. 

A complete implementation of all these effects does not exist to date. In our implementations 
l -^2 | 13 |^ have introduced a medium-induced term in the splitting function — item 1. above 
— and a formation time effect as an effective cut-off of the radiation — item number 4. above. 
The implementation of extra mechanisms is under way. 



2 Tiie modified splitting probability 

The central assumption in our approach is a factorization of the vacuum and medium contri- 
butions to the splitting probability P{z) as two independent and additive terms: 

Ptotiz) = Pvac(2) ^ i^tot(^) = P..ciz) + ^P{z,t,q,L,E). (1) 

This additivity have been proved at the level of the one-gluon inclusive spectrum, where the total 
contribution after computing the Feynmann diagrams of Fig. [T]can be factorized into a vacuum 
plus a medium term. By matching with the vacuum case, the corresponding modification to the 
splitting probability is written as 
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Figure 1: The three contributions to the gluon radiation in the presence of a medium (represented by a shaded 
box). The first two are included in the medium-induced gluon radiation term, while the third corresponds to the 

usual vacuum splitting. 



where dl'^'^'^ / dujdk\ is given by the medium- induced gluon radiation used previously for RHIC 
phenomenology ISI. This spectrum depends on two parameters, the medium length L and the 
transport coefficient q. The latter of these quantities encodes all possible information about the 
medium properties as temperature, density, etc. in a single parameter, with the interpretation 
of the average transverse momentum squared that the gluon gets per mean free path in the 
medium. This simple medium modification is implemented in the standard final-state showe ring 
routine PYSHOW in PYTHIA ^'^^ The corresponding routine is publicly available at the sit 

As an attempt to relate virtuality ordering in the PYTHIA parton shower and space-time 
development, we consider the formation time of the radiated gluons as tform = and 
subtract this quantity, after each splitting, to the total length of the traversed medium. This 
mixed approach allows us to implement space-time evolution and treat vacuum and medium 
radiation on the same footing through the modified splitting probabilities ([1]). 



3 Results 



The first expected effect of including an additive term to the splitting probability is an increase 
of the jet multiplicity. The typical medium radiation is also broadened in angle and softened in 
radiated energy with respect to the vacuum case. 

In Fig. [2] we show the default results at parton level, for two values of the transport coeffi- 
cient g = 5 and 50 GeV^/fm. We observe a suppression of high-z = log(l/2;)) particles and 
a big enhancement of particles with low-intermediate z-values, as expected. We also observe a 
suppression of high-pT particles and the corresponding enhancement of intermediate-py parti- 
cles. The pT-spectrum should be softer than vacuum at low transverse momentum since low-pr 
particles should be kicked towards higher values of the transverse momentum. However we find 
a clear enhancement of low-p^^ particles. Here, the lack of exchange of energy and momentum 
with the medium plus energy conservation in PYTHIA, is making a large effect. Finally, we see 
that the angular distribution broadens with increasing transport coefficient, as expected. 
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Figure 2: Intrajet parton distributions in ^ (left), pr (middle) and (right) for a gluon of initial energy Ejet = 100 
GeV in a medium of length L = 2 fm and for different transport coeffi cien ts q = (black), 5 (red) and 50 (blue 

lines) GeVVfm. Figure from Ref H^. 
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